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We study single-electron transport through a double quantum dot (DQD) monitored by a capaci- 
tively coupled quantum point-contact (QPC) electrometer. We derive the full counting statistics for 
the coupled DQD - QPC system and obtain the joint probability distribution of the charges trans- 
ferred through the DQD and the QPC consistent with the fluctuation theorem (FT). The system 
can be described by a master equation with tunneling rates depending of the counting fields and 
satisfying a generalized local detailed-balance relation. Furthermore, we derive universal relations 
between the non-linear corrections to the current and noise, which can be verified in experiment. 
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I. INTRODUCTION 

The well known fluctuation-dissipation relatior^i be- 
tween the conductance of a mesoscopic system and its 
noise is fundamentally important but limited to the low- 
bias, linear response regime. The fluctuation theorem 
(FT)^"— extends this relation to the high-bias and non- 
linear regime. It is formulated in terms of the entropy 
production in a mesoscopic system AS' during time r 
and relates the probabilities Pt{t^S) to flnd negative 
or positive values, T^-S*, in the simple relation 



-AS 



(1) 



The FT ID) has been verified in chemical, biophysi- 
caU^i^Sj and condensed matter— experiments. Re- 
cently it has also been tested in low-temperature elec- 
tron transport experiments with a serially coupled dou- 
ble quantum dot (DQD)--, and with an Aharonov-Bohm 
interferometer^-. In these cases the entropy production 
is related to the Joule heat AS = qeVs/T, where q is 
the number of electrons transmitted through the DQD 
during the measurement time r, Vs is the source-drain 
voltage bias. Hence the FT ([T]) takes the form 



Pr{q) 



exp 



'_qeVs\ 



(2) 



Experimental tests of the FT ([2]) in mesoscopic con- 
ductors are challenging because of the usually strong 
coupling to the environment, and, on top of that, the 
back action of the charge detector, which is relatively 
enhanced at low temperatures. In this paper we derive 
the full counting statistics (FCS) of the charge transport 
through a DQD capacitively coupled to a quantum point 
contact (QPC), which serves as a detector. We study the 
effect of both the electromagnetic environment and the 
back action of the QPC electrometer on the transport 
properties of the DQD. In earlier workS^ we had shown 



that the QPC back action leads to a violation of the FT 
^ when only the two-terminal DQD is considered. Here 
we explicitly show that the combined system of the DQD 
and the QPC satisfy a four terminal version of the FT 



Pr{q,q') 



= exp 



qeVs + g'eVqpc 



T 



(3) 



Here P{q,q') is the joint probability distribution for 
q charges to be transferred through the DQD and q' 
charges through the QPC during the time r, and Vqpc 
is the QPC bias voltage. We reveal that the FT dSj is 
the consequence of generalized local detailed balance re- 
lations, Eqs. (|28II30|) presented below, between the DQD 
transition rates, which turn out to depend on the count- 
ing field of the QPC. 

The direct experimental test of the relation is dif- 
ficult because the rate of electron tunneling through the 
QPC is too fast and the measurement of the joint proba- 
bility distribution Priq, q') is impossible. Therefore, sim- 
ilarly to Ref. [l^ we propose to test the relations between 
the non- linear transport coefficients, Eqs. psp and p9l) . 
These relations are the consequence of the FT ([3]) and 
can be verified by measuring the currents and the noises 
of the DQD and the QPC as functions of the bias voltages 
Vs and Vdqd- 

The FT has been previously discussed in Refs. 0-11 
at a formal level. It has been studied in detail for a sys- 
tem of two capacitively coupled quantum dot o^^i"'^'' , both 
of which are fully described by master equations. Here we 
explicitly derive it for the specific case of coupled DQD 
and QPC. Our analysis requires involved calculations an 
is based on the separation between the time-scales rele- 
vant for the dynamics of the DQD, which is a slow clas- 
sical system described by a master equation, and of the 
QPC, which is a fast quantum system, whose dynamics 
we describe in a path integral approach. As a result, the 
combined system is described by a master equation with 
tunneling rates modified by the QPC back action. We 
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take into account two back action mechanisms: (i) back 
action due to the switching of the QPC current between 
several values corresponding to different DQD charging 
states, and (ii) back action due to the electromagnetic 
environment excited by QPC shot noise. An approach 
similar to ours, which, however, ignores (ii), has been 
earlier developed by Braggio et al?^ 

The paper is organized as follows: In Sees. [TT] we 
present the experimental setup and our main results. In 
Sec, nil Al we specify the full microscopic model, then de- 
rive in Sec. HUB] the master equation for the DQD, and 
use it to find the PCS of the combined DQD + QPC 
system. Finally, in section HVl we summarize our results. 



(a-1) - 



II. MODEL AND MAIN RESULTS 
A. Setup 

We consider a DQD capacitively coupled to a QPC, 
which serves as a charge detector— i^l. The schematics 
of the system, as used in the experiment by Fujisawa et 
al^, is shown in Fig. [1] (a-1). Electrons tunnel between 
the two quantum dots and the leads via high-resistance 
tunnel junctions. The level spacing in the quantum dots 
is large enough such that only one level in each dot par- 
ticipates in the charge transport. The gate voltages ap- 
plied to the dots are tuned in such a way that only three 
charge states of the DQD contribute to the transport: 
the state with occupied left and empty right dots (state 

L)), the state with empty left and occupied right dot 
(state \R)) and, finally, the state with both dots occupied 
(state \D)), see Fig. [T] (a-2). We note that our analysis 
equally applies to the case when the states \L), \R), and 

0) (both dots empty) are occupied, i.e. in the vicinity of 
another tricritical point in Fig.. [T] (a-2). The state of the 
DQD is characterized by the vector p of the occupation 
probabilities of these states. 



= {pl,Pr,Pd) ■ 



(4) 



We assume spin degeneracy and from now on ignore the 
spin degrees of freedom. It is justified because the tunnel 
coupling between the leads and the dots is supposed to 
be small, and hence the non-trivial spin related physics, 
such as Kondo effect, is irrelevant. 

Due to the coupling to the DQD, the QPC transmis- 
sion and, as a consequence, the QPC current switches 
between three values corresponding to the DQD charge 
states |_R), \L} and \D), see Fig. [1] (b). We denote the 
corresponding QPC transmissions as T\r),T\l),T\d) and 
assume the following relation between them Tp^ < T\l) < 
7\u) ■ In addition, we assume that the QPC is tuned to 
the tunneling regime, where all transmission are small, 

T\n) < 1- 

Previous experiments revealed that the intrinsic capac- 
itive coupling between the QPC detector and the DQD 
leads to a strong parasitic back action of the former on 
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FIG. 1: (a-1) Setup used to count single-electron tunnel- 
ing events through the DQD in real time. The impedance 
^cxt(i^) models the electromagnetic environment which af- 
fects the QPC transport properties, (a-2) Low bias stability 
diagram of the double quantum dot. The solid lines indicate 
the equality of the energies of the neighboring charging states. 
The operating point is indicated by a black cross. It is cho- 
sen in the vicinity of a tricritical point, where the energies of 
the three charging states |L), |7i), and \D) are equal to each 
other, (b) The QPC current switches between three values 
corresponding to the three charge states of the DQD. (c) The 
relevant transition processes. Circles represent the double-dot 
states and arrows represent the directions of the transitions, 
(d) An equivalent circuit of the system. 



the DQD transport^"^'^*. The mechanism of this backac- 
tion is as follows: The fixed, i.e. time-independent, QPC 
bias voltage Vqpc is applied to and the QPC and the seri- 
ally connected environment impedance Z^^ti^) (see Fig. 
[T](a-1)). Since both of these elements produce noise, the 
instantaneous voltage drops across the QPC, WQPc(i)i 
and across the external impedance, Wext(Oj fiuctuate in 
time, while their sum stays constant, VQPc{t) -l-Wcxt(^) — 
Vqpc- Through the capacitive coupling, the fluctuations 
of the voltage across the QPC are transferred to the DQD 
and act in the same way as fluctuations of the gate po- 
tentials applied to each of the two quantum dots. These 
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fluctuations, in turn, modify the tunneling rates^S. 

The setup shown in Fig. [IJa-1) allows one to distin- 
guish the direction of the electron tunneling, and, hence, 
is suitable for testing the FT 1^ provided the back ac- 
tion from the QPC charge detector is minimized. For 
example, the switching of the QPC current between the 
levels \L) and \R) corresponds to the transfer of one elec- 
tron from the left to the right dot, while the switching 
\R) — >■ \L) signals a transfer in the opposite direction. 
Below we will show that the setup of Fig. [lJa-1) also 
allows testing the extended version of the FT 

B. Full counting statistics and generalized master 
equation 

In this subsection we summarize our main findings, 
while the details of the derivation can be found in Sec. lIIII 
The coupled system of the DQD and the QPC is de- 
scribed by the joint probability distribution Pr{q,q') of 
the charges q and q' being transfered through the DQD 
and the QPC, respectively. We introduce the correspond- 
ing characteristic function (CF) 

Zr{x,x)^Y.Pr{q,q')e'^'^''''', (5) 

I 



and the cumulant generating function (CGF) 



T{X,x)= lim -lnZ,(A,x), (6) 

r— foo 7" 

where A and x are the counting fields of the DQD and 
QPC respectively. 

The charge transport though the DQD occurs by 
sequential tunneling of electrons through the high- 
resistance tunnel junctions. This process is described by 
a simple master equation for the occupation probabilities 
of the three DQD charge states. Since we are interested 
in the full counting statistics of the electron transport, 
we need to keep track of the counting fields A and x- In 
this case the occupation probabilities Q depend on the 
counting fields (see Eq. ([5^ for precise definition), and 
the master equation acquires the form 

5,p(r)-M(A,x)p(T), (7) 

where M is a 3 x 3 transition matrix. 



M 



'TRLiOj-TDLiO) 

TB.L{x)e'^ 

^DLix) 



TLBix)e-'^ 
-rLfl(O) - Tdb.{0) 
^dr{x) 



^LDix) 

^rd{x) 

-rLD{o)-rRD{o) 



Tl + (5 J"o 
+ (5 Jo 

I'd + (5 Jo 



.(8) 



From the solution of the rate equation ([7]) the CF ([5]) is 
found as follows 



n=L,B.,D 



(9) 



The outlined procedure differs from the usual one^^ only 
by the dependence of the tunnel rates r„in on the QPC 
counting filed x- 

Let us now discuss various parameters which appear 
in the transition matrix ([8|). The parameters J„ are the 
CGF of the QPC at a fixed charging state of the DQD, 
indicated by the subscript n. They are given by 



Tu 



\n) 



eVqpc 



e'x - 1 



(10) 



The correction (5Jo(x) is defined in Eq. (175)) and is due 
to influence of the environment on the QPC. 

The transition rates in the matrix ([5]) are defined as 



follows 

rLfl(x) 
^rl{x) 

^dr{x) 



RD 



DL 



(x) 



(x) = 



2TTtlPdd^^{AR ~ Al) , (11) 
2TTtlPM^^{AL - Ar) , (12) 

Tl J dEf{-AR - + E)Pl,^{E), (13) 

Tl J dE[l - fi~AR ~^iL~ E)]Pl,^{E), 

(14) 

dEf{-AL - fiR + E)Pr^^{E), (15) 



Tld{x) ^ dE[l~fi^AL~m~E)]PR..xiE)- 

(16) 

Here we introduced the hopping matrix elements be- 
tween the two quantum dots tc] the energy differences 
A„ = En — Ee) between the energies En [n — L, R, D) of 
the DQD charging states; the tunnel coupling strength 
between the dots and the leads Tlj^r', and the Fermi 
function f{E) = l/(e^/^ + 1). The back action of the 
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QPC on the transport properties of the DQD is encoded 
in the functions P„^(£^), where n — dd,L,R indicates 
one of the three tunnel barriers in the DQD. These func- 
tions describe both the effect of the QPC shot noise and 
the excitation of the QPC environment^^'^° and depend 
on the QPC counting field x- They are defined as 



dt 
2^ 



(17) 



The dimensionless parameters = —kr = i^ddl"^ char- 
acterize the capacitive coupling between the DQD and 
the QPC and are found from the equivalent circuit shown 
in Fig. [11(d): 



Cl{Co + Cg + Cr) 



+ {Co + Cl){Cg + Cr) + 2CoC 



(18) 



where C = 2Cg + Cl + Cr. The function reads 



(19) 



$x(i) = ie^ J dLu\Zt{uj)\^ 



where 



UJ" 

- (5o(^) + e2-^o(x)) (5o(-c^) + e2-^o(x))] • (20) 

Here ^o(x) is defined by Eq. ([TU]) with 7[„) replaced by 
the QPC transmission 7" in the regime when both quan- 
tum dots are empty, and 

Zt{uj) - [-iuC -f e^TA + l/^cxt(w)] (21) 

is the effective impedance seen by the QPC with C being 
the total effective capacitance seen by the QPC (see Fig.[T] 
(d)), 



C 



CL[3Cg + Cg{Cg + Cr) + 2Co(2Cg + Cr)] 
iCl + {Cg + Cl){Cg + Cr) + 2CoC 



(22) 



S^{uj) is the non-symmetrized current noise spectral 
power given by the sum of the QPC shot noise and ther- 
mal noise of the environment, >S'x(w) = Si^^i'-^) + 'S'em('^), 



5/,x(w) 



TT 



{uj - eV^QPc)e 'X (tj + eVQPc)e^^ 



1-e- 



1 - e 



5*0111 (t^) = Re 



^ext(w) 



2w 



(23) 
(24) 



C. Generalized detailed balance and fluctuation 
theorem 

The tunneling rates (|llj|16p obey the generalized de- 
tailed balance relations. In order to derive them, we first 



note that in equilibrium the noise of the QPC environ- 
ment satisfies the usual detailed balance relation 



S'om(w)e 



-;9w 



5*0 



(25) 



where /3 = 1/T is the inverse temperature. Since the 
QPC noise depends on the counting field x it obeys a 
generalized detailed balance relation of the same type. 



(u;)e-''- = Si,^{-u,). 



(26) 



With the aid of these two relations one can prove the 
identity for the modified Pn.xiE) functions: 



Pn-x+il3eVQrciE)e 



-I3E _ 



(27) 



These identities combined with the properties of the 
Fermi function lead to the generalized local detailed bal- 
ance relations between the tunneling rates: 

Trl{x) - riflX-X + «/3e"t^QPc)e''(^^-^«\ (28) 
Tdl[x) ^TLD{-X + ifieVQvc)e^^'"'^'"'\ (29) 
Tdr{x) = r;j,D(-x + z/3eFQPc)e^('^«+^-), (30) 

The relations (I28ll30p constitute the first main result of 
our paper. They differ from the usual local detailed bal- 
ance relations in one important point: in our model the 
tunneling rates depend on the counting field x which has 
to be transformed as x — ^ —X + */3eVQpc to ensure the 
detailed balance symmetry. 

It is known^ that the local detailed balance is sufficient 
to prove the FT for an isolated DQD system We will 
now demonstrate that the relations (|28II30I) lead to the 
FT © for the coupled system of DQD and QPC. First, 
we note that the CGF of the QPC (fTO)) has the property 



^n{x) = ^n{-X + «/3eVQPc), 



(31) 



which is equivalent to the FT for an isolated QPC ([2|). 
The same applies to the correction induced by the envi- 
ronment ((75)) . 



5Mx) - SM-x + ifieVQvc). 



(32) 



Next, in the long-time limit r 3> the character- 

istic function ([S|) approaches the asymptotic exponential 
behavior 



Z.(A,x)«e-^(^''^) 



(33) 



where the CGF J' is the eigenvalue of the matrix ([5]) with 
the largest real part. It is to be found from the equation 



det [TE - M{X,x)] =0. 



(34) 



With the aid of the identities (1281132^ one can prove the 
following property of the determinant 

det [TE - M(-A + i(3eVs, ~x + tPeVQPc)] 

= det[J"£;-M(A,x)] ■ (35) 



5 



Hence the function x) obeys the symmetry 

^(A, x) = -F(-A + il3eVs, -x + ^PeV^Pc). (36) 

Since in the long-time limit the characteristic function 
P3p possesses the same symmetry, one can easily ver- 
ify that the joint probability distribution, given by the 
inverse Fourier transformation 



Pr{q,q') 



dXdx 



~i\q-ixq' 



Zr{X,x), (37) 



satisfies the FT in the form ([3]). The proof of the FT is 
the second important result of our paper. 



D. Experimental verification of tiie FT for the 
coupled system 



DQD charge states A^. These parameters may be in- 
dependently measured, which provides additional check 
both for the FT and for our model. The corresponding 
expressions read 



T 



5^/ _ a^DQD 



dvl 



dVs 



g3^g/3(A^+A„)|p 



LRX ld^rd 



e RD 

(AB)2 

2 ^LD dvs^RD + ^RD dy^TLD 
LR 



^LR^LD^RD 

2Tr 



where we have introduced the parameters 



(41) 



Here we are going to discuss how one can verify Eq. Q 
in the experiment. The direct measurement of the joint 
probability distribution Pr{q,q') is almost impossible 
with current experimental techniques because the tunnel- 
ing of electrons through the QPC is a very fast process. 
Therefore, we propose a different route which requires 
only current and noise measurements. Namely, we pro- 
pose to check the universal relations between nonlinear 
transport coefficients imposed by the FT^iSsAS. The first 
one relates the lowest-order nonlinear conductance of the 
DQD with the correction to its noise and reads 



T 



dvl 



Vs=Vqpc=0 



dVs 



(38) 



^^s=Vqpc=0 



Here / is the current through the DQD, and iSdqd is 
the DQD current noise. The identity (|38p can be derived 
from both the two terminal ([2]) and the four terminal ^ 
versions of the FT—. 

One can derive a relation analogous to Eq. psp for the 
QPC current and noise. However, since the nonlinearity 
of the QPC I — V curve is weak and easily smeared by 
noise and temperature, the latter relation is more difficult 
to verify in the experiment. 

The second relation, which we propose to test, reads 



2T 



dVsdVQPc 



as. 



QPC 



^^s=Vqpc=0 



dVs 

dSnl 



Vs=Vqpc=0 



i9Vqpc 



(39) 



Vs=yQPc=o 



Here 



5„i = / dt{{l{t) - (/))(/QPc(0) - (/qpc))) (40) 



is the non-local cross correlation of the QPC and DQD 
noises. 

Our model allows us to express all the derivatives, 
which appear in the Eqs. P8I39|) in terms of the tun- 
neling rates Pmn and the energy differences between the 



B 



e^^^TLDTRD + TLRiTLD + Trd). (42) 



The nonlinear response of the QPC current reads, 



2T 



dVsdVc^PC dVs 

TT 

Tlr{Tld - Trd) + TldTrd{2 + e'^^n) 



A^B 



^ ^j, T\R) -7[j) ^g(2A^+AH) 
TT 

Tlr{Tld - Trd) - TLDTRDef'^^il + 2e-^^^) 



A^B 



(43) 



In Eqs. (j41l43l) all tunnel rates are taken at zero QPC 
counting field, x = and zero bias voltages Vqpc — Vs — 
0, and can be directly measured in experiment. Finally, 
we find that the cross correlation of the noises equals to 
zero. 



dS^\/dVQPc\y^^ 



Vs=Vqpc=0 



= 0. 



(44) 



This result significantly simplifies the test of the relation 
P9p making rather complicated measurement of the non- 
local noise unnecessary. 



III. DERIVATION OF THE MASTER 
EQUATION 

In this section we provide a detailed derivation of the 
master equation ((7|) and the expressions for the tunneling 
rates ()imi6|) . For simplicity, here we put the electron 
charge equal to 1, e = 1, and ignore electron spins. 
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A. Model Hamiltonian 

We describe the system shown in Fig. [1] (a-1) by the 
Hamihonian 

H = H,,s + Hi2 + Ht, (45) 

H,es = Hl+Hr + Hi+H2+ iJcm + -H"dQD • (46) 

Here the Hamiltonians 

Hr = ^ Crfc CLrk'^ark , {r = L,R) , (47) 
k 

describe the electrons in the left and right leads coupled 
to the DQD, while the Hamiltonians 

Hj=J2^jka,k^a,k, (j = l,2), (48) 

k 

describe the leads 1 and 2 on both sides of the QPC. The 
Hamiltonian of the double quantum dot Hdqd reads 

Hbqd = J2 En\n}{nl (49) 

n=L,R,D 

where En is the energy of the DQD charge state \n). 

The electron transport through the QPC is described 
by the tunnel Hamiltonian 

Hi2 = u ^ (ti2e'^a2fe^aife. + /i.e.) , (50) 

kk' 

where ti2 is the hopping matrix element, 
u= ^ Un\n){n\ 

n=L,R,D 

is the operator which accounts for the effect of the DQD 
charge state on the QPC tunneling, and 9 is the phase in- 
duced by voltage fluctuations in the electromagnetic envi- 
ronment. If both quantum dots are unoccupied the QPC 
transmission probability reads T — 47r^i^2'^i^^2j where 
Vi,V2 are the densities of states in the leads 1 and 2. If 
either left or right quantum dot is occupied the tunnel 
matrix element is suppressed and the transmission prob- 
ability changes to T\n) = u^T, n = L, R, D. 

The tunneling between the quantum dots and the left 
and right leads is described by the tunnel Hamiltonian 

r=LM k 

-Hice"^""^ + /I.C., (51) 

where f = L/R for r — R/L. Due to the capacitive cou- 
pling between QPC and DQD the phase 9 also appears in 
this Hamiltonian. Finally, the electromagnetic environ- 
ment with the impedance Zcxt(w) is modeled by a set of 
LC resonators in the spirit of Caldera-Leggett model^. 
It is described by the Hamiltonian 

a 



where the charge operator Q is related to the phase op- 
erator by a commutation relation [0, Q] = j; ip^ are the 
coordinates, or phases, of the environment LC-oscillators 
and the corresponding charges are defined by the com- 
mutators [(pa,Qa] — i- The spectral density of LC- 
resonators is chosen to reproduce the impedance of the 
electric circuit Zcxt(w) 

Y] i^i^ - ^a) - + ^o)\ -Re , y (53) 



B. Full counting statistics 

In this subsection we derive a formal expression for the 
CF ^ applying the FCS approach for the reduced den- 
sity matrix"^^. We split the system into six reservoirs: five 
electronic ones, namely the leads 1,2,L,R and DQD, and 
the bosonic reservoir of LC-oscillators. The reservoirs are 
described by the Hamiltonians iJi, H^, Hl^ Hr, ^^dqd 
and Hem respectively. At negative times i < all reser- 
voirs are decoupled and the system is described by the 
factorized density matrix 

Pa = PL Pk Pi P2 Pcm J2 PniO)\n){n\, (54) 

n=L,FI.,D 
g~(ff,-M.Afr)/T 

P"- ^ xr (e-(^--''.^-)/^) ' ^^^^ 

where Nr — '^/^ark^o.rk {r = 1,2) are the parti- 

cle numbers in the r-th reservoir, fj,r are the chemical 
potentials with the chemical potential of the bosonic en- 
vironments being equal to zero, /^om = 0, and p„ are the 
initial occupation probabilities of the DQD states \n). 

At time t — the tunnel Hamiltonians Hi2,Ht are 
turned on and at i = r they are turned off again. We 
introduce the initial nI^\n[^^ and final n'i^\n[^^ num- 
bers of the electrons in the leads L and 1, and define the 
probability Pr{q,q') for q electrons to tunnel out of the 
reservoir L and q' electrons - out of the lead 1 as follows 

Pr{q,q') - J2\(f\^''"^\')\'('\Po\^) 

if 

where \i) and |/) are eigen states of the hamiltonian H^cs- 
The CF dni) is then expressed as, 

if 

= Tr(e-'-f^+>oe'-^-^) . (57) 

Here we defined 
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The only difference between these Hamiltonians and the 
original one (|45p is the additional phase shift in the tun- 
neling matrix amplitude. Namely, after the transforma- 
tion ([55]) the tunnel Hamiltonians (|50l5ip acquire the 
form 



Hi 



2± 



kk' 



+ /I.C.) , (59) 



H' 



T± 



{tL e^'^ e'''-''^^aLk^\R){D\ + h.c. 



aRk'\L){D\+h.c. 



(60) 



Defining the Keldysh contour K , which starts at i = 0, 
propagates to t = r and then returns back to i = (sec 
Fig. [2]), one can interpret as the system Hamiltonian 
on the forward branch of this contour, while — its 
Hamiltonian on the backward branch. 

Next, we split the Hamiltonians ((58|) into two parts 
H± — Hq± + Ht±, where the "non-perturbed" Hamilto- 
nian consists of the reservoir part and the QPC hopping 
part, 



Ho 



± 



Hr, 



Hi 



2± ■ 



Switching to the interaction representation we define 
HT±{t)i = eiHo±t ^Q-iHo±t rewrite CF in the 
form 



Zr = Tr 



TkgxpJ^-z jdtHT±{t)^ e^^°-^e-*^«+>o 



(61) 



where Tk is the time ordering operator along the Keldysh 
contour. By formally introducing the identity operator 
Sn=L D I"") ("I ~ 1' expressed as the sum over the 
charge states of the DQD, into this formula we transform 
Zr to the form Eq. ([SJ with the occupation probabilities 



Vn{\x,T) = Tr 



TKe^v[-i JdtHT±{t)i 



xe'""-^\n){n\e-'"«+''po 



(62) 



which satisfy the normalization ^^PniO, 0; t) — 1. 

We find the probabilities ^ performing the expansion 
in powers of the tunnel Hamiltonian Ht± and using the 
real-time diagrammatic technique^"—. The details of 
this procedure are outlined in the next subsections. 



C. Zeroth order. 

If Ht± — the probabilities p„ read 

pi°nT) = 7r„(T,0)p„(0), (63) 



H+ 



H- 



t = 



FIG. 2: Keldysh contour K. The time runs from the right 
to the left. The total system evolves with the Hamiltonian 
H+ for the forward branch (upper contour) and for the 
backward branch (lower contour). 



where 



7r„(T,0) =Tri 



Pi P2 Pc 



(64) 



is the system propagator at fixed DQD charge \n). The 
Hamiltonians Hn± = Hi + H2 + i^em + ("■|-ffi2± 1*^) in- 
clude the capacitive coupling between the DQD and the 
QPC. Employing the standard technique'^^ we write the 
propagator 7r"(r, 0) in the form of the path integral over 
the fluctuating phase 9 



7r„(T,0) = jv[e] exp(i5QPc[0, 



Un]+iS,M , (65) 



The action of the QPC is derived by tracing out elec- 
tron degrees of freedom in leads 1 and 2 and subsequent 
expansion of the result in powers of the tunnel matrix ele- 
ment ti2 (see Appendix |A|). The result of this procedure 
is the nonlinear action, which makes the path integral 
([65)1 impossible to evaluate exactly. To overcome this ob- 
stacle we assume that the impedance of the electromag- 
netic environment is sufficiently small, Zs{uj) <C e^/27r, 
at all relevant frequencies. In this regime the action 
Sqpc[(),u] becomes quadratic in phase fluctuations and 
takes the form, 



'QPC 



[0, 



dtdt' Si^^{t,t')u{t)u{t') 



K 



It carries the information about the DQD charge state |n) 
through the function u{t), which, in general, depends on 
time and switches between the discreet values m„. The 
function Si^^{t,t') is related to the QPC current noise 
and has the 2x2 matrix structure in Keldysh space 



(67) 



where the superscript + indicates the forward branch of 
the Keldysh contour , while — stands for the backward 
one. In the limit of long measurement time r — >■ 00 
Si^x{t,t') depends on the time difference t — t' , which 
allows one to make the Fourier transformation. After- 
wards the greater and lesser components of this matrix 
take the form are expressed with Eq. 



(68) 
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Within our approximation the causal and anti-causal 
components do not depend on the counting field 



±p. 



27ri uj ~ uj' 



(69) 



where the second part is the Cauchy principle value. The 
action of the electromagnetic environment reads 

Son.m= I dt^e^+'- I dtdt'd{t)S,^{t,t')9{t'), (70) 

where we have already traced out Q, Qa and (fa- 
The lesser/greater component of the matrix Scm equals 
S^^iu:) = S'em(T'^) [see Eq. ([M]) ]. while causal/anti- 
causal components are given like Eq. ([S^ . 

Combining all these components, we express the total 
action ^eff = S'qpc + S'cm in the form 

5eff[0,u] - 5qpc[o,u] + ]- I dtdt'e{t)D-\t - t')e{t'), 



(71) 



where 



C 



-u? 

^2 



5++(c.) -5+-(c.) 
-^-+(c.) 5--(c.) 



«-^o(x) 



1 
1 



(72) 



Jx)(x) is given by Eq. (fTO|) with Tn replaced by T, and 

S~)^ = "5*/,^ + 5*0111 • ('''3) 

Note that we have split the action (|71l) into the contribu- 
tion, which depends only the DQD charging state, i.e. on 
M(t), and the contribution, which depends only the fiuc- 
tuating phase 0, but its dependence on u(t) is ignored and 
we put u=\. This approximation is justified as long the 
modulation of the QPC transmission probability by the 
DQD is smaU \Tn-T\<-T. 

4 components of the full photon matrix Green function 
are, 



|Zt(^)p5T±(c.) 



(74) 



(75) 



where the function ^^(iS) and the impedance Zt{uj) are 
defined in Eqs. (PH]) and (PT|) respectively. 

Let us now evaluate the first term of the effective action 
([7T|). In the limit r -4 cx) we get 



«'S'qpc[0, Ur, 



2 s+-{0) + sy+io) ~ 5++(o) - sr-{o) 



(76) 



which is the CGF of the QPC Hence the propa- 

gator (p5|) takes the form 

7r„(T,0) w / P[6']e^^'=« = e^i^r.{x)+STo(x)]^ (77) 



where 6To{x) oc — lndet[D -^j is the correction to the 
QPC cumulant generating function, which reads 



(78) 



The correction ST{x) is caused by weak the Coulomb 
blockade effects in the QPC shunted by the impedance 
Zcxti^)^^- By virtue of Eq. ((26)) and a similar identity 
S'Z^e"^" = S'+r one can show that fi, 



Hence SJ'{x) satisfies the FT 



X+«/3Vqpc — 



D. Second order: Photon assisted tunneling 

Let us now discuss the second order expansion in Hj^ . 
Since the Gaussian effective action ([7T|) is equivalent to 
some effective bosonic bath, we adopt the corresponding 
diagrammatic rules proposed in Ref. 3^. Figure [3] shows 
few typical diagrams. In these diagrams the time runs 
from the right to the left; the horizontal lines on the up- 
per and lower branches of the Keldysh contour denote 
the states of the DQD; the directed solid lines, connect- 
ing the two branches, depict the Green functions of the 
reservoirs. The Fourier transformed lesser and greater 
components of the latter Green functions read 



j^^{Lo)=rrf{T^±tir), {r = L,R). 



(79) 



where the rates F^ — 2'Ktr^Vr determine the coupling 
between the dot r and the lead r. 

The sum of the two diagrams (a-1) and (a-2) reads 



(a) _ 



dt2dhty''^ (r, t2,tu0) e'(^«-^-)(*^-*i)pL(0), 



(80) 



where the propagator is defined as follows 



(a) _ { AHl^Ii -iKadS^-iHn-iti-T) -iHRj^(T-t2) 



V[e] exp l^S^s \ 0, u'^"^]^ + ijds J(°) (s) 0{s) 

:^ g''SQPc[o,«<°']+r5J^o(x) 

xexpj^-^ J dsds' J^''\s)D^{s,s')J^''\s') 



The source J^°^ and the potential u^"'' are expressed via 
the delta function 6 and the step function 9 defined on 
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(a-l) |R> 



(a-2) |R) 



>f — -■ 

dd 'l 
'2 ,-' 


m 




*" ^ * 

dd 

~ \. « 


m 


= T t = 

(b-1) |R) |D> 


(b-2) |R> |D) 










15(0) 


|D) |R> 


(c-2) |D> |R> 




m 






m 


(d+) |L) |R) 


(d-) |R> 


dd 


m 




dd 

— » x' »• i< 





FIG. 3; Typical diagrams contributing second order expan- 
sion in terms of the tunnel couplings, x represents hopping 
between the two dots tc- Horizontal lines represent the dot 
state on the Keldysh upper and lower branches. Directed solid 
lines are reservoir lines. Dotted lines are correlation functions 
of the Gaussian phase fluctuations denoted by e"^'''**'^ etc. in 
the text. 



the Keldysh contour K: 

J^'Hs) = Kdd{S{s,t2+)-S{s,h^)), (82) 
u'^"^\s) =. UL + {uR-UL){e{ti-,s)^e{sM+))- (83) 

Provided ^ T, eV , we may use Markov approxima- 
tion and put 



iSi 



QPC 



:Fr{x) {T-t2)+FL{x)ti- 
The remaining double integral reduces to 

dsds' J^''\s)D^{s,s')J^''Hs') 



(84) 



K 



D--(h,ti) + D++{t2M) 



ii^d^xi^i -*2), 



(85) 



where the function $^(0 defined in Eq. (|19p . We 
note that at x = the function $o(0 coincides with the 
phase correlation function which appears in the so called 
P(£;)-theory (see e.g. Eq. (18) of Ref. IH). Physically 
this function describes the back action of the QPC charge 
detector on the DQD through the photon assisted tun- 
neling^: caused by the QPC shot noise. 

Thus after adopting Markov approximation the dia- 
grams (a-l) and (a-2) give the result 

5p^^\t) « fdt e-*^"(x)+^«(x) (r-t) Tj,i^{x)pf\t) ,(86) 



Likewise the diagrams (b-1) and (b-2) give the contri- 
bution 



(87) 



where the factor e containing the DQD counting field 
comes from the lines describing the electron tunneling 
through the DQD. 

The diagrams (c-1) and (c-2) describe the opposite 
transition processes and their contribution can be ob- 
tained from that of the diagrams (b-1) and (b-2) by in- 
terchanging the DQD states D ^ R and changing the di- 
rection of the reservoir line corresponding to 7^^(aj) 
7^~(— cj). In this way one arrives at the expression (fT3)) 
for the tunneling rate F^i^. 

Let us now demonstrate that the function Pn.x{E) sat- 
isfies the generalized detailed balance ([?7|) . It follows 
from the following property of the phase correlation func- 
tion 'I'x(^)- 



*-X+i/3VQPc(i) 



duj- 



S'o(w) - e"''^*5'_x-i-i/3VQpc 
duj- 



X (5o(c^)+^o(x)-e'"^*+^^^^xH 



(88) 



where we used, |Zt(w)| = \Zt{—uj)\ andJ7^(cj) = r2^(— w) 
in addition to Eqs. (|26|) and ([GS)) . The detailed balance 
relation (|27p can now be verified. 

Let us now consider the diagrams (d-|-) and (d-), which 
define the diagonal tunneling rate. The contribution of 
these diagrams can be evaluated as follows 



-tt 



dU 







t2 



dtl^(''+)(T,<2,<l,0) 
di2 7r('^-'(T,t2,tl,0) 



dti 
^0 

«(iJ^-iJH)(t2-*l)p^(0), 



(89) 

where for the propagator 7r^''^\ are given by path inte- 
grals with the Gaussian action Soft with the source fields 

J^'^^Hs) = km{5{sM±) - 5{sM±)) , (90) 
u('^±)(s) =UR+ [UL - ub) i0{h±,s) - 0(s,t2±)).(91) 

Thus the sum of the diagrams (d-|-) and (d-) is evaluated 
to the expression 



where TrlIx) is defined by Eq. ([T 



(92) 
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where we neglected x dependence of Tlr since it gives 
only small correction to J'r{x)- 

Other rates T^n are found analogously and also obey 
the generalized local detailed balance relations. We note 
that at X = the master equation ^ reduces to the 
ordinary one, which conserves the probability p„ = 1. 
Higher order expansions are systematically calculated in 
the same way. In this way, we obtain the integral form 
of the master equation Q. 



IV. SUMMARY 



In conclusion, we have derived the full counting statis- 
tics for the coupled system of a double quantum dot 
(DQD) and a quantum point contact (QPC). The non- 
trivial part of our analysis is the development of the the- 
oretical framework which accounts for both the quantum 
charge transport through the QPC and the classical se- 
quential tunneling through the DQD. We find that the 
combined system should be described by a generalized 
master equation, in which the tunneling rates are mod- 
ified by the QPC back action and depend on the QPC 
counting field. We prove the generalized detailed bal- 
ance relations between the tunneling rates corresponding 
to time-reversed processes. These relations, in turn, lead 
to the fluctuation theorem for the joint probability dis- 
tribution of the charges transfered through the DQD and 
the QPC. Finally, we derive universal relations between 
the non-linear corrections to the currents and the noises, 
which follow from the fluctuation theorem and can be 
verified in the experiment. 

We would like to emphasize that in this paper we have 
limited ourselves by the regime of weak tunneling in the 
DQD. At stronger coupling between the quantum dots 
and the leads higher-order quantum corrections have to 
be taken into account and the master equation ([7]) is not 
valid any more. In this case one has to deal with the FCS 
of non-equilibrium strongly correlated systems (see e.g. 
Ref. [srt ). and we leave that problem for the future. 

We thank Toshimasa Fujisawa, Hisao Hayakawa, Mat- 
tias Hettler, Kensuke Kobayashi, Bruno Kiing and Keiji 
Saito for helpful discussions. This work has been sup- 
ported by Strategic International Cooperative Program 
of the Japan Science and Technology Agency (JST) 
and by the German Science Foundation (DFG) and the 
Okasan-Katoh foundation. 



Appendix A: Effective QPC action 

The Keldysh action describing the QPC in Eq. (p5|) is 
defined as follows 

•^^ [j=1.2 k 

-ti2u{t)^ ^a2k{ty aik' (t) e' 



te{t)+x{t)/2 



+C.C.]}, 



(Al) 



Here the QPC counting field xi^) equals x on the for- 
ward, and —X on the backward branches of the Keldysh 
contour. Integrating out the Grassmann fields ajk and 
taking care of the proper boundary condition^, we ar- 
rive at the result 



'QPC 



= -iTvG 



XI 



(A2) 

^-^'''^ [-h2it,t') g-\t,t') )■ ^^^^ 

Here the times t and t' run over the whole Keldysh con- 
tour K, and the trace is understood as the time integral 
along K. The components of the matrix (jASp read 



h2{t,t') = ti2ii(i)e*«(*)+^'^«/2^(i,t'), (A4) 
gj{t,t') = ^g,fe(t,0, (A5) 



(A6) 



Expanding the action (|A2[) up to second order in ti2, we 
obtain the Ambegaokar-Eckern-Schon action^^ 



QPC = « dtdt' u{t)u{t')tl29i{t,t')g2{t',t) 

JK 

^^i[e(t)-0(t')]+i[x{t)-x(t')]/2 _ 



(A7) 



Further expanding this action in 9, we get the action (j66p 
with 

Si,x^{t,t') = u{t)u{t')tl^gi{t,t')g2{t',t) 

xe*[xW-x(«')]/2 ^ ^ ^^3^ 

Appendix B: Nonlinear transport coefficients 



In order to derive Eqs. (j4T|) and (|43|) we solve the 
equation for the CGF ([M)) perturbatively in the count- 
ing fields X, A and bias voltages Vs,Vqpc, and use the 
definitions 



>DQD 



5, 



QPC 



9^ 



5, 



nl 



dXdx 
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Several additional symmetries make the calculations sim- 
pler. First, at Vqpc = 0, the local detailed balance for 

0) are 



the DQD tunneling rates (Eqs. (|28II30|) with x 
satisfied. Therefore at Vs — Vqpc = the following 
relations hold 



Tom 

dVs 
dVl 



dVs 



f^^^LpiO) 



2/3 



dVs 



dVs dVs 



dVs 



LP 



(0) 



Second, as long as we keep only the corrections oc 
the dependence on x is determined by the CGFs for the 
QPC (fTU]) . which in this limit reads 



Tn « 7T„)/(27r) [FQPc(^x) + Tiixf 



(Bl) 



The dependence of the tunneling rates on x can be ne- 
glected, = $0 since at low bias, the QPC contribution 
is dominant. 

Derivatives of the transition rates over Vqpc show how 
strongly the tunneling in the DQD is affected by the QPC 

= 



shot noise. One can check that dvQpcSi.o{(^) 
and thus the first derivatives vanish, i.e. 



Vqpc=0 



(0) 



Vqpc=0 



= 0. 



(B2) 



By utilizing these results, we arrive at Eqs. (HTI) and P5)) 
after expanding the characteristic equation (j34p in pow- 
ers of Xj Vs and Vqpc and comparing the coefficients 
(A more systematic derivation for the higher cumulants 
in Ref. is also applicable). 
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